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ABSTRACT 

Context. Intermediate-mass stars are an important ingredient of our Galaxy and key to understand how high- and low-mass stars form 
in clusters. One of the closest known young intermediate-mass protoclusters is OMC-2 FIR 4, located at a distance of 420 pc, in 
Orion. This region is one of the few where the complete 500-2000 GHz spectrum has been observed with the heterodyne spectrometer 
HIFI on board the Herschel satellite, and unbiased spectral surveys at 0.8, 1, 2 and 3 mm have been obtained with the JCMT and 
IRAM 30-m telescopes. 

Aims. We aim to disentangle the core multiplicity and to investigate the morphology of this region in order to study the formation of 
a low- and intermediate-mass protostar cluster, and aid the interpretation of the single-dish line profiles already in our hands. 
Methods. We used the IRAM Plateau de Bure Interferometer to image OMC-2 FIR 4 in the 2-mm continuum emission, as well as 
in DC0+(2-l), DCN(2~I), C'^S(3-2), and several CH3OH lines. In addition, we analysed observations of the NH3(1,1) and (2,2) 
inversion transitions made with the Very Large Array of the NRAO. The resulting maps have an angular resolution which allows us 
to resolve structures of 5", equivalent to ~ 2000 AU. 

Results. Our observations reveal three spatially resolved sources within OMC-2 FIR 4, of one or several solar masses each, with hints 
of further unresolved substructure within them. Two of these sources have elongated shapes and are associated with dust continuum 
emission peaks, thus likely containing at least one molecular core each. One of them also displays radio continuum emission, which 
may be attributed to a young B3-B4 star that dominates the overall luminosity output of the region. The third source identified 
displays a DCO*(2-l) emission peak, and weak dust continuum emission. Its higher abundance of DCO^ relative to the other two 
regions suggests a lower temperature and therefore its possible association with either a younger low-mass protostar or a starless core. 
It may alternatively be part of the colder envelope of OMC-2 FIR 4. 

Conclusions. Our interferometric observations evidence the complexity of the intermediate-mass protocluster OMC-2 FIR 4, where 
multiple cores, chemical differentiation, and an ionised region all coexist within an area of only 10000 AU. 

Key words. ISM: individual objects: OMC-2 FIR 4 - ISM: molecules 

1 . Introduction In the past decades, most of the attention has been devoted 

to the study of low-mass and high-mass star formation, reveal- 
Star formation starts in dense fragments within clouds of molec- ing differences between the two which a re not only dynamical 
ular gas and dust and ends when nuclear hydrogen burning be- (Beuther et al. 2007 Zin necker & Yorke H^T), but a lso chem- 
gins at the centre of a newly formed star. During this evolu- jcal in nature (Bottinelli l2006l Ceccarelli et aL l2007l l. The for- 
tion, the physical processes and different phases involved are mation of intermediate-mass (IM) stars represents the bridge be- 
not necessarily the same for stars of all masses. Indeed, in the tween low-mass and high-mass star formation and can bring 
field of star formation there has traditionally been a division important pieces of information to understand how differently 
between low-mass stars (M < 2 Mq), intermediate-mass stars the process of star formation works in the two mass extremes. 
(2 Mo < M < 8 Mo), and high-mass stars (M > 8 Mq). As their high-mass counterparts, IM stars emit a non negligible 

amount of radiation in the far ultraviolet (FUV), and are thus 

* Based on observations carried out with the IRAM Plateau de Bure ^^le to photodissociate and ionise their suiTOundings. Therefore, 

Interferometer. IRAM is supported by INSU/CNRS (France), MPG the study of IM star formation can aid the cuiTent debate on 
(Germany) and IGN (Spain). 
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the formation mechanism of high-mass stars (see Zinnecker & 
Yorke l2007l for a review). 

While low-mass stars may form in isolation or in loose stel- 
lar aggregates, high-mass stars are always seen to be formed in 
densely clustered environments. IM star formation also occurs 
in clustered mode (e.g. Fuente et al. l2007l l and has been found to 
be the transition between loose low-mass aggregations and tight 
high-mass clusters (Testi et al. ll999l l. IM protoclusters represent 
excellent laboratories to study clustered star formation at almost 
the full range of stellar masses. Observationally, they can typi- 
cally be found closer {d < I kpc) than high-mass protoclusters, 
and they are also less complex. 

Despite all this, relatively little work has focused on the early 
stages of IM star formation. Only a few individual sources have 
been studied in detail (e.g. Fuente at al.'2009' Grimier et al. 120091 
I20T0I Sanc hez-Monge et al. 2010 Palau et al. 201 1 , van Kempen 
et al. l2012l l. and therefore much remains unknown about the for- 
mation and first evolutionary stages of these objects. 

The present work focuses on the IM protocluster OMC- 
2 FIR 4, with a luminosity that has been reported to be between 
50 L0 (Adams et al. 2012 ) and 1000 Lq (Grimier et al. 2009) . It 
is located in the Orion A complex, 2 pc North of the famous 
Orion Nebula (M42) and the Trapezium OB association, and 
lies ~420 pc away from the sun (Menten et al. 120071 Kim et 
al. 12008). 

Due to its relative proximity, multiple studies have been car- 
ried out towards this region and the OMC-2/3 filament where 
it is located. These include continuum observations at several 
wavelengths (e.g. Chini et al. 1997 Johnstone & Bally 1999, 
Nielbock et al. 120031 ). and various molecular line studies (e.g. 
Tatematsu lT993l Cesar oni & Wil sonlT994l Taka hashi e t al. l2508l 
Tatematsu et al. l2008l Liu et al. ISOTTI Li et al. l20T2l ). Most of 
them, however, are coarse angular-resolution observations, very 
useful to investigate the overall structure of the whole filament, 
but not sufficient to resolve the substructure within the OMC- 
2 FIR 4 region. The only dedicated interferometric observations 
towards this source are Very Large Array (VLA) maps carried 
out at 3.6 cm by Reipurth et al. (1999), and a recent study con- 
ducted with the Nobeyama Millimeter Array (NMA) at 3 mm by 
Shimajiri et al. (120081 hereafter S08). The latter reveal several 
cores and the authors speculate that an external outflow driven 
by the nearby FIR 3 is responsible for triggering fragmentation 
in FIR 4, highlighting the complexity of this region. 

OMC-2 FIR 4 is also one of the targets of the HerscheQ 
Guaranteed Time Key Programme Chemical HErschel Surveys 
of Star forming regions (CHESSq Ceccarelli et al. 2010 Kama 
et al. l2012l ). which has conducted unbiased spectral surveys with 
the HIFI spectrometer (de Graauw et al. 120101 ) of eight star 
forming regions, each representative of a particular aspect of the 
star formation process: evolutionary stage, mass of the forming 
star, and/or interaction with the surroundings. The HIFI half- 
power beam width ranges between 11" at 1900 GHz and 41" 
at 500 GHz, thus covering most, if not all, the area of OMG- 
2 FIR 4. High-angular resolution observations are therefore a 
crucial tool to help disentangle the information hidden in the 
HIFI line profiles. 

With this in mind, we report in this paper 2-mm contin- 
uum and line observations performed with the IRAM Plateau de 



Bure Interferometer (PdBI) and ammonia maps obtained with 
the VLA towards OMG-2 FIR 4. We describe these observations 
in Sect. |2j in Sect. |3]we present the maps and spectra obtained, 
which we discuss and interpret in Sect. |4] Finally, Sect. |5]sum- 
marises the main conclusions of this work. 

2. Observations 

2.1. Plateau de Bure observations 

We observed OMG-2 FIR 4 with the IRAM Plateau de Bure 
Interferometer (PdBI) in two tracks on 15 October 2010 and 15 
April 2011, using the most compact D and G configurations of 
the array. The data presented in this paper cover a spectral win- 
dow of 1.8 GHz obtained with the WIDEX correlator, whose 
channel spacing is 1.95 MHz (~4 km s '), centred at a fre- 
quency of 143.4 GHz. The phase centre of the observations is 
Q'(J2000) = 05*'35'"26.97P, 5(J2000) = -05°09'56.77". The sys- 
temic velocity of OMG-2 FIR 4 is Vlsr = 1 1 .4 km s^' . 

The bandpass of the receivers was calibrated by observing 
0923-H392 in Oct 2010 and NRAO150 in Apr 2011, with a 
flux density of 3.55 Jy and 2.53 Jy, respectively. Amplitude and 
phase calibrations were achieved by monitoring 0430-1-052 and 
0528-1-134, whose flux densities were determined to be 1.32 and 
1.67 Jy, respectively, at 143.4 GHz. The uncertainty in the am- 
plitude calibration is estimated to be 30%. 

The data were calibrated and analysed with the GILDASJJ 
software package developed at IRAM and the Observatoire de 
Grenoble. We obtained continuum maps from the line free chan- 
nels and subtracted the continuum from the line emission di- 
rectly in the (u,v)-domain. The resulting 1 cr RMS and beam size 
for the naturally weighted continuum map are 4.9 mJy beam"' 
and 4'.'87 x 2'.'73, respectively (see TablelTJ. The maps presented 
in this paper are not corrected for the primary beam, but we did 
apply the correction to measure the reported fluxes. 

2.2. Very Large Array ammonia observations 

The Very Large Array (VLAQ) was used to observe the am- 
monia {J,K) = (1,1) and (2,2) inversion transitions simultane- 
ously (project AG556). Twenty adjacent pointing centres, cov- 
ering OMG-2 and OMG-3 regions, were carried out on 29 July 
2000 and 24 September 2000, with the array in the D (compact) 
configuration. The FWHM of the primary beam at the frequency 
of the observations is approximately 2'; the field centres were 
separated by 1' to ensure full sampling. Each field was observed 
for a total on-source integration time of approximately 30 min- 
utes. The absolute flux scale was set by observing the quasar 
1328+307 (3G286), for which we adopted a flux of 2.5 Jy. The 
quasars 0605-085 and 0539-057, with bootstrapped fluxes of 
2.4 Jy and 1.0 Jy, respectively, were observed regularly to cali- 
brate the gains and phases. Bandpass calibration was performed 
by observing the bright quasar 0316-1-413 (3G84) with a flux 
density of 14.2 Jy (for July) and 17.5 Jy (for September). 

The data reduction followed the VLA standard guidelines 
for calibration of high-frequency data, using the NRAO package 
AIPS. The NH3 (1,1) and NH3 (2,2) fines were observed simul- 
taneously in the four IF correlator modes of the VLA (with two 



' Herschel is a European Space Agency (ESA) space observatory 
witli science instruments provided by European-led principal inves- 
tigator consortia and with important participation from the National 
Aeronautics and Space Administration (NASA). 

^ http://www-laog.obs.ujf-grenoble.fr/heberges/chess/ 



3 The GILDAS package 

http://www.iram.fr/rRAMFR/GILDAS 



available 



at 



* The Very Large Array (VLA) is operated by the National Radio 
Astronomy Observatory (NRAO), a facility of the National Science 
Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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Table 1. Parameters of our continuum observations and those of Shimajiri et al. (I2008I I 



Parameter 



Tliis work 



Sliimajirietal. (20081 



Interferometer 
Frequency (GHz) 
Weigiiting 
Primary Beam (") 
Baseline range (k/1) 
Beam size (arcsec) 
PA. of tliebeam (°) 
R.M.S. (mJybeam'') 



PdBI 

143.4 
natural 

35 

10-80 

4.87x2.73 

-0.76 

4.9 



PdBI 

143.4 
uniform 

35 

10-80 

3.74x2.74 

-12.98 

4.6 



NMA 

92 
natural 

82 

2.9- 115 

6.96x4.03 

-34.59 

1.2 



NMA 

92 

uniform 

82 

2.9-115 

6.55x3.34 

^0.31 

1.4 



polarisations for each line), providing 63 channels with a spac- 
ing of 0.62 km s ' across a bandwidth of 3.125 MHz, plus an 
additional continuum channel containing the central 75 per cent 
of the total bandwidth. The bandwidth was centred at the sys- 
temic velocity Vlsr =11.4 km s " ' . 

Calibrated uv-data were loaded in the software MIRIAD 
(Sault et al. 1995) to perform the imaging. While the continuum 
emission was found to be negligible (we measured an upper limit 
of about 2 mJy at the frequency of the observations) we sub- 
tracted any possible related contribution from line-free channels 
using the task UVLIN. We combined and imaged the data from 
the central part of the mosaic, covering the OMC-2 region, with 
the task MOSSDI. The rms noise in the naturally weighted maps 
is 5 mJy beam"' per 0.62 km s"' channel. We applied a Mv-taper 
function of 60 kA (3" in the image plane) to improve the signal- 
to-noise ratio. The resulting synthesised beam is 6'.'2 x 5'.'3, with 
PA=-5°. 



3. Results 

3. 1 . Millimetre continuum emission 

Figure [T] shows our PdBI naturally and uniformly weighted 2- 
mm continuum maps (white contours). For comparison, these 
are overlaid on the 3-mm continuum images obtained by SOS 
with the Nobeyama Millimeter Array (NMA). From their uni- 
formly weighted image, SOS identify 1 1 cores at 3 mm within 
OMC-2 FIR 4, whose positions according to their Table 4 are 
marked with black crosses in the figure. The parameters of our 
continuum observations and those of SOS are listed in Table [T] 

The PdBI images allow us to resolve two elongated struc- 
tures in OMC-2 FIR 4, separated by ~ 5", surrounded by more 
extended emission which covers a total size of about 15". The 
uniform map suggests the presence of additional fragmentation 
which cannot be completely resolved at the angular resolution 
of the observations. On the other hand, the NMA images dis- 
play fainter, more complex and more extended emission, which 
is largely below the 3cr level (first contour) in the case of the 
uniform map. 

It is worth mentioning that the primary beam size of our 
PdBI observations is 35", considerably smaller than that of the 
NMA maps, which is S2". Similarly, the shortest baseline at- 
tained with the NMA is 2.9 kA while that of the PdBI is 10 kA. 
This results in a higher sensitivity of the NMA to extended struc- 
tures. More specifically, our shortest baseline implies that our 
observations are sensitive to structures smaller than 17" down to 
a level of 1 0% (se e Eq. AS in Wilner & Welch^li or Eq. A6 in 
Palau et al. 120101) . while for the NMA images this number is as 
high as 56". Thus, we are not sensitive to the cores identified by 
SOS that fall outside the outermost PdBI contour In spite of this, 
we still have several concerns regarding the core identification 
made by SOS. First, it was based on the uniform map, which has a 
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Fig. 1. Interferometric continuum maps of OMC-2 FIR 4 with 
natural (top) and uniform weighting (bottom). White contours 
represent the 2-mm continuum emission obtained with the PdBI, 
while black contours, as well as the background colour scale, 
correspond to the 3-mm continuum emission seen with the NMA 
(SOS). Contours for both the PdBI and NMA maps start at 3cr 
and increase by steps of 3cr. The bottom-left ellipses repre- 
sent the respective beam sizes. The white open circles mark 
the primary beam FWHM for the PdBI observations. The black 
crosses mark the NMA 3-mm cores as reported in Table 4 of 
SOS. 



low signal-to-noise (S/N) ratio; second, the coordinates given for 
some of the cores do not coincide with the corresponding emis- 
sion peak positions in the image; and third, two pairs of cores lie 
within a beam size of each other We believe the NMA natural 
map is more reliable than the uniform one given its higher sen- 
sitivity. This already shows a number of clearly separated cores, 
some of which appear to be followed also by the PdBI contours. 
In conclusion, even though we prefer to be cautious about the 
reality of the 1 1 cores reported in SOS, we are convinced that 
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Table 2. Lines observed with the PdBI and VLA interferometers 



Line 



v/' 
(GHz) 



(K) 



(s-') 



Itr RMS" 
(mjy beam"') 



CH,0H(3i 3 - 


- 2,,2 A+) 


143.865801 


28.3 


1.07x10-' 


CH3 0H(3o,3 - 


- 2o,2 E) 


145.093707 


27.1 


1.23x10-5 


CH30H(3-i 3 


- 2-1,2 E) 


145.097370 


19.5 


1.10x10-5 


CH3 0H(3o,3 - 


- 2o,2 A+) 


145.103152 


13.9 


1.23x10-5 


CH30H(32,2 - 


- 22,1 A-) 


145.124410 


51.6 


6.89x10-'^ 


CH3 0H(32,i - 


- 22,0 E) 


145.126190 


36.2 


6.77x10-'^ 


CH30H(3-2 2 


- 2-2,1 E) 


145.126392 


39.8 


6.85x10-'^ 


CH30H(3i,2 - 


- 2,,i E) 


145.131855 


35.0 


1.12x10-5 


CH30H(32,i - 


- 22,0 A+) 


145.133460 


51.6 


6.89x10-'^ 



3.7 



DCO+(2-l)^ 


144.077289 


10.4 


2.12x10-' 


3.6 


e'*S(3 - 2) 


144.617101 


13.9 


5.78x10-' 


3.2 


DCN(2 - 1)^ 


144.828002 


10.4 


1.27x10-* 


3.4 


NH3(1,1) 

NH3(2,2) 


23.69450 
23.72263 


23.3 
64.4 


1.68x10-' 
2.24x10-' 


5.0 
4.7 



'"' Millimetre line frequencies from the CDMS catalogue: http://www. astro.uni-koeln.de/cdms/] 

'** — indicates no map is presented in this paper 

'^' Transition with 6 hyperfine components 



OMC-2 FIR 4 is composed of several cores. Based on our PdBI 
maps, we will hereafter refer to the dominant eastern region as 
main source, and to the western one as west source. 

The continuum fluxes, measured inside the 5cr contour of the 
corresponding maps, are 1.06 + 0.32 Jy and 1.10 + 0.33 Jy for 
the natural and uniform maps, respectively. We can estimate the 
mass of the emitting clump by considering that all the flux is 
emitted by dust. This is a reasonable assumption, since Reipurth 
et al. ( [19^ reported a 3.6-cm flux of 0.64 mJy for OMC-2 FIR 4 
(see Sect. 14.31 1. and therefore the contribution of free-free emis- 
sion at 2 mm is negligible, i.e. about three orders of magnitude 
below the total flux we measure. Assuming that the dust emis- 
sion is optically thin, the mass, Mdust is given by 



Mdu 



Sy(f 



K, B,{Ti) Ri 



(1) 



where 5 y is the flux we meaure inside the 5cr contour at 2 mm, d 
the distance to OMC-2 FIR 4, T^ the dust temperature (adopted 
to be 50 K from SOS), ByiTi) is the Planck black body func- 
tion for a temperature Td, and R^ is the dust-to-gas ratio, as- 
sumed equal to 0.01. /Cy = K(){—i^ is the frequency-dependent 
dust mass opacity coefficient. Chini et al. d 1 9971) derived /3 - 2 
from the Spectral Energy Distribution fit towards OMC-2 FIR 1 
and FIR 2. Since FIR 4 is in the same molecular filament as 
FIR 1 and FIR 2, we have adopted/? = 2, and kq - I cm^ g-' at 
vo = 250 GHz (Ossenkopf & HenningfT994b. 

The resulting mass, from the naturally-weighted 2-mm con- 
tinuum map, is 9 + 3 Mq. This should be considered as a lower 
limit to the total mass of OMC-2 FIR 4 (dense cores plus enve- 
lope), since a percentage of the total flux is expected to be filtered 
out by the interferometer. We can estimate the fraction of flux 
lost with the aid of the work published by Crimier et al. (120091 ). 
who modelled the Spectral Energy Distribution (SED) function 
observed towards OMC-2 FIR 4 at infrared and sub-millimetre 
wavelengths. The flux they predict at 2 mm is 2 Jy, implying a 
flux loss of 45% in our PdBI continuum observations. 

Crimier et al. (12009) derived a mass of 30 M© for the enve- 
lope of OMC-2 FIR 4, i.e. higher than our estimate by a factor 3 . 
While this difference may be partly due to the loss of flux in our 
interferometric maps, we note that a different choice of dust tem- 
perature in Eq.[T]can significantly change the obtained mass. For 



example, adopting Td = 20 K, which is a reasonable dust tem- 
perature for the source outer envelope (Lis et al. 1998. J0rgensen 
et al. i2006. Crimier et al. 12009) , yields a mass of 26 Mq. 

3.2. Line emission 

The list of detected transitions are presented in Table |2] together 
with the Icr RMS values of the coiTesponding velocity channel 
maps (see Sects. [mi and l3T2!2b . Out of the nine methanol lines 
detected, eight are blended due to the poor spectral resolution 
of the observations. Therefore, the methanol maps shown in this 
paper correspond to the only isolated transition we observe, i.e. 
CH3OH(3,l,+0 - 2,l,+0) at 143.87 GHz, whose upper level en- 
ergy is 28 K and which will hereafter be refeiTed to as CH3OH 
for simplicity. 
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Fig. 2. PdBI spectrum obtained towards the phase centre of 
OMC-2 FIR 4 with the WIDEX correlator. The species of the 
detected lines are labelled. 



3.2.1. PdBI molecular line spectra and maps 

The 1.8-GHz wide spectrum obtained at the peak coordinates of 
OMC-2 FIR 4 with the WIDEX correlator at PdBI is presented 
in Fig. |2] where the detected species are labelled. Since we are 
also in possession of IRAM 30-m single-pointing data towards 
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Fig. 3. Velocity-integrated PdBI contour maps obtained towards OMC-2 FIR 4, overlaid on the naturally weighted continuum map 
(grey scale). The bottom-left ellipse in each panel represents the beam size. Left: CH30H(3, 1,-1-0 - 2,l,+0) (black) and C^'*S(3-2) 
(magenta). Centre: DCN(2-1) (black) and DCO^(2-l) (magenta). AU the contours start at 3cr and increase by steps of 3cr. The Icr 
RMS values are 0.16, 0.14, 0.17, and 0.08 Jy beam' km s', respectively for the CH3OH, C^^SO-l), DCN(2-1), and DCO+(2-l) 
integrated maps. Right: Polygons representing the three differentiated components, main, west, and south (see text), used to extract 
average spectra. The white crosses mark the positions of the three sources identified in this work (see Table |4). 



the same source (Lopez-Sepulcre et al., in prep.), we provide 
in Table [3] an estimate of the velocity-integrated flux loss for 
each PdBI line. This ranges from 13% for CH3OH to 73% for 
DCO^(2-l), and indicates the presence of an important contri- 
bution from extended emission in these tracers which is filtered 
out by the interferometer 

Table 3. Flux loss for the PdBI molecular fines 



Line 


Flux loss 




(%) 


CH3OH 


13 


DCO+(2-l) 


73 


€^"8(3 - 2) 


60 


DCN(2 - 1) 


72 



Given the low spectral resolution of the WIDEX spectrum 
(1.95 MHz, i.e. 4 km s '), it is not possible to analyse the shape 
of the line profiles. However, this is beyond the scope of the 
present work, which aims primarily at investigating the spa- 
tial distribution of the emission in the diff'erent tracers. To this 
end we provide, for each millimetre molecular tracer, velocity- 
integrated maps in Fig. [3] and velocity channel maps in Fig. |4] 
The synthesised beam size of these images is 4.8x2.7 arcsec, and 
is depicted on the bottom-left corner of each panel in the figures. 

If we focus first on the integrated maps (Fig. |3]l, the morpho- 
logical complexity of the source becomes evident from a com- 
parison of the emission distribution in the different tracers: 

- The CH3OH emission is dominated by the main source (see 
Sect. 13.11 above), which displays an elongated shape, and 
shows hints of a secondary component towards the west (the 
above mentioned west source), mimicking the spatial dis- 
tribution of the 2-mm continuum emission. This morphol- 
ogy does not change when seen in the blended CH3OH lines 
around 145 GHz. 

- C''^S(3-2) displays a more compact emission that traces both 
the main and west sources. 

- The DCN(2-1) emission peaks at the main source, but it ex- 
pands towards the west and the north, following the elon- 
gated morphology of the west source as seen in the contin- 
uum images. 



- DCO^(2-l) peaks to the South of the main and west sources, 
a different region which we name south source, not well 
traced by any other line. 

In summary, from the integrated maps we conclude that 
each molecular line traces the region diff'erently. We diff'eren- 
tiate three spatially separated emission components, which we 
have named main, west, and south sources, lying ~ 0.01 pc (i.e. 
2000 AU) apait from each other These regions are schematically 
depicted in the right panel of Fig. |3]by three magenta polygons 
that have been manually defined to pinpoint similar central ar- 
eas for each while keeping approximately their morphological 
shapes. From these we have extracted the respective continuum 
and line fluxes for each tracer The resulting fractional fluxes 
with respect to the total flux, measured inside the 3cr map con- 
tour for the lines and the 5cr contour for the continuum map, are 
listed in Table |4] The central coordinates of each source corre- 
spond to the emission peaks of the 2-mm continuum image for 
the main and west sources, and those of the DCO^(2-l) map 
for the south one. Notice that, with the exception of DC0^(2- 
1), due to the relatively small area of each polygon, the sum 
of fluxes is smaller than the total flux in the whole region. The 
sum amounts to 136% in the case of DCO^(2-l) partly because 
the polygons corresponding to the main and south region over- 
lap around the position of the peak emission, and partly because 
there is an amount of contributing flux outside the 3cr contour. 
This table reflects the predominance of the main source's flux 
over that of the other two for all the tracers except DCO^(2-l), 
for which the south source dominates the overall emission. 

Using Eq. [1] we have computed the mass, M, of each source 
for two values of the dust temperature: 20 and 50 K. These 
are listed in Table |5] together with the corresponding source - 
averaged gas column densities derived from the equation: 



A^H 



M 



II mu 



Q.(f 



(2) 



where m^ is the mass of the hydrogen atom, /i = 2.33 is the 
mean molecular mass in units of hydrogen atom masses, and Q 
is the solid angle subtended by the dust continuum emission. 
Assuming the depth of the region equals the projected one (i.e. 
5000 AU), we derive an average density in the range 0.9 - 2.6 x 
10^ cm-3 for OMC-2 FIR 4. 
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Table 4. OMC-2 FIR 4 sources: coordinates and fractional fluxes for each tracer 



Source 


R.A.(J2000) 


Dec.(J2000) 


CH3OH 


C^''S(3-2) 


DCN(2~ 


-1) 


DC0+(2- 


-1) 


Cont 


(natural) 


Total" 






13.0 


3.7 


6.3 




1.2 






1.1 


Main(%) 
West(%) 
South(%) 


05''35"26.97^ 
05''35™26.63^ 
05''35"'26.80» 


-05°09'57.8" 
-05°09'56.2" 
-05°09'59.3" 


49 
15 
13 


51 

25 
11 


61 
30 
12 




39 
32 
65 






35 
18 
12 



'°* Measured within the 5cr and 3cr contour for the continuum and line maps, respectively; flux units in Jy for continuum and Jy km s ' for lines 
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Fig. 4. Channel contour maps of each line detected with PdBI towards OMC-2 FIR 4. Contours start at 5cr and increase by steps of 
5cr. The Icr RMS values are listed in Table |2] The velocity of each channel is marked on the top-left corner of each panel. Black 
crosses mark the phase centre of the observations, i.e. a(J2000) = 05^35™26.971\ 5(J2000) = -05°09'56.77", while red crosses 
mark the positions of the three sources identified in this work (see Table lU. The bottom-left ellipses in each panel represent the 
beam sizes. 



The channel maps in Fig. |4] reveal more details about the locity channels (i.e. covering a width of about 18 km s '), and 

spatial distribution of the molecular line emission. Despite the therefore some kinematical information can already be extracted 

low spectral resolution of the observations, the emission in all from these maps. We note, however, that in the case of DCN(2- 

the tracers but DCO^(2-l) is detected above 5cr across 5 ve- 1) this velocity spread is partly due to the hyperfine splitting 
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(see Sect. 13.2.31 below). In Fig. |5]we show the average spectra 
within the main, west, and south sources, defined by the poly- 
gons shown in Fig. [3] (right panel). Figure |6] presents the first 
moment maps for the four molecular tracers observed with the 
PdBI (colour scale), with overlaying contours of the correspond- 
ing velocity-integrated maps. The spectra, channel maps and ve- 
locity maps reflect the following properties : 

Table 5. OMC-2 FIR 4 sources: masses and column densities" 



1 - 



Source 


M(20 K) 

(Mo) 


M(50 K) 

(Mo) 


NH2i20 K) 
(10" cm-2) 


(10 


,(50K) 
23 cm-2) 


Total 


25.7 


9.2 


19.2 




6.9 


Main 
West 
South 


9.1 

4.7 
3.0 


3.2 
1.7 
1.1 


29.0 
18.7 
124 




10.4 
6.7 

4.4 



'"* Derived from the naturally-weighted 2-mm continuum map 



The main source displays a velocity gradient in the south- 
west to north-east direction. 

The west source as seen by CH3OH and C^'*S(3-2) appears 
slightly blue-shifted with respect to the systemic velocity of 
the source. 

The DCO^(2-l) emission is significantly the narrowest in 
line width and concentrates almost exclusively on the south 
source. 

A secondary, well-separated condensation is distinguished in 
the CH3OH and DCN(2-1) maps to the south-west of the 
main source, suggesting the presence of an additional phys- 
ical source. The DCO^ emission also extends towards this 



While the channel maps suggest the presence of more 
sources apart from the three we have identified so far, we cannot 
confirm it with 100% certainty at this spatial resolution, and even 
less guess how many more are present. We therefore conclude, as 
stated above, that there are at least three distinct sources within 
OMC-2 FIR 4, leaving open the possibility of further multiplic- 
ity within each of these three regions. 

3.2.2. Ammonia emission 

The NH3(1,1) and (2,2) line spectra seen towards the peak 
NH3(1,1) emission coordinates are shown in Fig.|7] with the dif- 
ferent hyperfine components marked in red. The limited band 
width of the correlator allows us to cover only one pair of satel- 
lites for each transition, apart form the main component. In the 
case of NH3(2,2), only the main component is detected. 

Velocity-integrated ammonia maps are shown in Fig. |8] 
where the emission corresponds to the main hyperfine compo- 
nent of the NH3(1,1) and NH3(2,2) transitions. These are over- 
laid on the 2-mm continuum PdBI image (grey scale) for com- 
parison. We stress that the presence of extended emission which 
was filtered out by the VLA results in a relatively high 1 cr RMS 
level in the clean maps. Indeed, Li et al. (120121) recently pub- 
lished these ammonia data merged with single-dish maps ob- 
tained with the Green Bank Telescope (GBT), and a comparison 
with our VLA-only spectra indicates a flux loss of about 40%. 
Nevertheless, in terms of morphology, the NH3 emission is well 
defined around OMC-2 FIR 4 and peaks to the South of the main 
source. The emission appears extended, with a roundish shape, 
and covers an area which is larger than that of the 2-mm contin- 
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Fig. 5. Spectra in the three main sources identified in OMC- 
2 FIR 4 (see text) for each line tracer. The vertical scale is the 
brightness temperature of the line, Tl, averaged across the area 
of the corresponding polygon (see Fig. |3). 



For completeness, in Fig.|9]we present the channel maps for 
the two NH3 transitions, with the velocity of the channels given 
on the top-right corner of each panel in km s ' . Unlike the mil- 
limetre molecular spectra (Sect. 13.2. U . the NH3 lines are rather 
narrow, as can be seen from the velocity range covered by the 
maps (see also the spectra in Fig. |7]and the FWHM values in 
Table[8]l. This may partly be due to the better spectral resolution 
and lower sensitivity of these observations. Notice indeed that 
the channel spacing of the VLA observations (~0.6 km s ') is 
substantially smaller than that of the WIDEX PdBI spectra. 

The channel maps show extended emission around the cen- 
tral coordinates of OMC-2 FIR 4, and several other sources of 
emission: two towards the noith-east and the south-west, where 
FIR 3 and FIR 5 are located, respectively (green circles), and a 
third one towards the south-east at the systemic velocity of the 
source (11. 4 km s"'), also seen in the single-dish 1.3-mm maps 
carried out by Chini et al. (II997) . 

3.2.3. Derivation of pliysical parameters 

Both the DCN(2-1) and DCO+(2-l) lines are split into six hy- 
perfine components. For DCO^(2-l), all of them lie within a 
velocity range of 0.3 km s"', while for DCN(2-1) they span 8 
kms"'. As a result, at the spectral resolution of our observations, 
they are all blended and unresolved, and therefore a hyperfine 
analysis is not possible. We have thus performed a Gaussian fit 
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Fig. 6. Velocity maps (colour scale, in km s ') of each line detected with PdBI towards OMC-2 FIR 4, with overlaying contours 
of the corresponding velocity-integrated maps. Contour levels are as in Fig. |3] The bottom-left ellipse in each panel represents the 
beam size. Crosses mark the phase centre of the observations, i.e. q'(J2000) = 05^35™26.97P, 5(J2000) = -05°09'56.77". 
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Fig. 7. VLA ammonia spectra obtained towards the peak coor- 
dinates of OMC-2 FIR 4, with a channel width of 0.62 km s '. 
The different hypeifine components are marked with red vertical 
lines. 

to each line and adopted a range of Tex between 10 and 100 K to 
estimate a range of column densities from: 



A^mnl = 



Q 



jT^dv 



c' Aui gue^ -I J(T..) - J(ng) 



(3) 



where v, c, h, k, A^u Tpeak, AV, Q, g^, E^p and J Ti^dv are, re- 
spectively, the frequency of the transition, the speed of light, 
Planck constant, Boltzmann constant, the spontaneous emission 
rate, the optical depth at the line intensity peak, the line FWHM, 
the partition function, the upper state degeneracy, the upper level 
energy and the velocity-integrated line intensity. This equation 
assumes optically thin conditions, and J(T) - -p^jwri, with 
Jbg = 2.7 K. We have applied the same approach for the C^'*S(3- 
2) transition. 

The results of the Gaussian fits and the derived molecular 
column densities for each source are shown in Table |6] where it 
can be seen that we are indeed dealing with optically thin con- 
ditions. Some of the rows in this Table probably coiTespond to 
"artificial" sources, since not all the sources we have identified 
are well-defined or traced by all the lines, as explained above. 
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Fig. 8. Velocity-integrated VLA ammonia contour maps ob- 
tained towards OMC-2 FIR 4, overlaid on the naturally weighted 
2mm-continuum map (grey scale). Contours for all the maps 
start at 3cr and increase by steps of 3cr. The Icr RMS values are 
20 mJy beam-i km s^'for both NH3(1,1) and NH3(2,2). White 
crosses mark the positions of the three sources identified in this 
work (see Table|4]l. 



We nevertheless present them for comparison and completeness. 
The abundances of each molecule in the different sources are 
also listed. They have been obtained from the ratio of A^moi to the 
gas column density, A^. For the latter, we have used the values 
Usted in Table|5] adopting both Td = 20 K and T^ = 50 K. 

As for the methanol lines, we analysed them by means of 
a Large Velocity Gradient (LVG) code, originally developed by 
Ceccarelli et al. (2003) and modified to include the methanol 
collisional coefficients with H2 by Pottage et al. (12004b . For this 
analysis, we used the methanol lines obtained by averaging the 
emission over the entire OMC-2 FIR 4 clump. The line intensi- 
ties thus measured are listed in Table |7] where the values shown 
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Table 6. OMC-2 FIR 4 sources: Gaussian fit results, column densities and abundances for each PdBI tracer 



Tracer 


Source 


(K) 


Vlsr 
(kms"') 


FWHM 
(kms-') 


Tex" 

(K) 




r 


A'mol 

(10'- cm-2) 


Abu. (20 K)* 

(10-'-) 


Abu. (50 K)* 
(10-'-) 


C3''S(3-2) 


Main 
West 
South 


0.37 (0.01) 
0.38 (0.01) 
0.15(0.01) 


10.9(0.1) 
9.3(0.1) 
9.3 (0.2) 


11.3(0.2) 
10.5 (0.2) 
10.7 (0.6) 


10 - 100 
10 - 100 
10 - 100 


0.055 
0.040 
0.017 


- 0.004 

- 0.003 

- 0.001 


15.3-38.2 
8.8 - 22.2 
4.0-10.3 


5.3-13.2 
4.7-11.9 
3.2-8.3 


14.7 - 36.7 
13.1-33.1 
9.1-23.4 


DCN(2-1) 


Main 
West 
South 


0.55 (0.01) 
0.67 (0.01) 
0.29 (0.01) 


10.6(0.1) 
9.8(0.1) 
9.0(0.1) 


10.0(0.1) 
9.8 (0.2) 
10.2 (0.4) 


10 - 100 
10 - 100 
10 - 100 


0.085 
0.052 
0.022 


- 0.006 

- 0.003 

- 0.001 


6.8 - 22.0 
3.7-12.5 
1.8-6.0 


2.3 - 7.6 
2.0-6.7 
1.5-4.8 


6.5-21.2 
5.5- 18.7 
4.1- 13.6 


DC0+(2-l) 


Main 
West 
South 


0.19(0.01) 
0.20 (0.01) 
0.38 (0.01) 


10.8 (0.2) 
11.1 (0.2) 
10.3(0.1) 


6.1 (0.3) 
6.3 (0.3) 
6.3 (0.2) 


10 - 100 
10 - 100 
10 - 100 


0.024 
0.020 
0.044 


- 0.002 

- 0.001 

- 0.003 


0.6 - 2.3 
0.6-2.1 
1.3-4.5 


0.2 - 0.8 
0.3- 1.1 
1.0-3.6 


0.6 - 2.2 
0.9-3.1 
3.0- 10.2 



*°' Adopted values 
**' Abundance determined using the dust temperature indicated in braclcets to derive the total gas column density (see Table [5) 
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Fig. 9. VLA NH3(1,1) and NH3(2,2) channel maps towards OMC-2 FIR 4. Contours start at 5cr and increase by steps of 5cr. The 
Icr RMS values are 5.0 and 4.7 mJy beam-' for NH3(1,1) and NH3(2,2), respectively. The velocity of each channel is marked on 
the top-right corner of each panel in units of km s' . Red crosses mark the positions of the three sources identified in this work (see 
Table|4]i. Green circles mark the position of FIR 3 (NE) and FIR 5 (SW) according to the 1.2 mm observations of Chini et al. (I1997I I. 
The bottom-right ellipses in each panel represent the beam sizes. 



for the 145-GHz lines do not correspond to single transitions but 
to ensembles of blended methanol lines around the indicated fre- 
quency and transition. We ran a large grid of models to cover a 
wide parameter space: N(CH30H) from lO'-' to 4 x 10'^ cm-^, 
density from 10"* to 10^ cm--' and temperature from 10 to 200 K. 
The observed lines were compared with the LVG model predic- 
tions and the best fit solution was found minimising the x~ value 
with respect to the three above parameters, assuming a source 
size of 15" X 15", as measured in our maps for the whole OMC- 
2 FIR 4 clump. We note that for the blended lines the predicted 
intensities were obtained by adding the intensity by each single 
line in the velocity interval. 

We found that the CH3OH column density is very well con- 
strained between 0.6 and 1 xlO'^ cm"^, roughly consistent with 
the value of ~ 2 x lO'"^ cm-^ reported by Kama et al. (1201 Oi l 
from the early results of the HIFI data. The density is between 
0.2 and 1 xlO' cm--', a range of values that matches the average 
density we obtain from the millimetre continuum emission (see 



Sect. 13.2. 11 1. The predicted optical depths of the CH3OH lines 
are low, around a few times 10-^. 

The observed lines do not allow to constrain the tempera- 
ture, except that it is larger than 20 K. The larger the CH3OH 
column density the larger the gas temperature. For example, if 
N(CH30H) is 8 x lO''* cm-^ then the gas temperature is larger 
than 50 K; for N(CH30H) = 1 x 10'^ cm-^ the gas temperature 
is larger than 120 K. More methanol transitions will be neces- 
sary to better constrain the temperature, but our analysis already 
shows that the gas emitting the methanol lines must be relatively 
dense. The best fit solutions of the LVG modelling have a re- 
duced ;t'^ below 0.5 and are summarised in Table |7] 

Finally, we performed a hyperfine line fit to the two ammo- 
nia spectra in Fig.Qusing the standard NH3(1,1) and NH3(2,2) 
methods in CLASS. We subsequently derived the coiTesponding 
rotational and kinetic temperatures, as well as the NH3 column 
density, following the steps summarised in Busquet et al. (20@ 
and assuming a beam filling factor equal to 1 . The results are re- 
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Table 7. CH3OH line intensities and results of the LVG analysis 



CH3OH line intensities 


Transition 


V 

(GHz) 


J TssdV 
(Kkms-') 


Notes 


3i,3 - 2i,2 A+^ 
3-1,3 - 2-1,2 E+ 

3-2,2 - 2-2,1 E+ 


143.866 
145.097 
145.126 


5.2+1.6 

28.8 + 8.7 

13.9 + 4.2 


isolated 
3 blended lines 
5 blended lines 


Results of the LVG analysis 


N(CH30H) 

(10'"* cm-2) 


Density 

(10* cm"') 


Temperature 
(K) 




6 
8 

10 


>2 

>2 

2- 10 


20-50 
>50 
> 120 





ported in Table|8] where the peak synthesised beam temperature, 
TsB, the velocity, Vlsr, and the optical depth, T^ain refer to the 
main component of the transition. We note that the errors on the 
temperature and column density estimates are of the order of or 
greater than the values themselves, mostly due to the large eiTors 
on TsB and r. Using the H2 column densities listed in Table |5] 
we obtain an NH3 abundance slightly above 10 '°. While this 
value is about 2 orders of magnitude smaller than what is typi- 
cally reported in the literature, there are studies of ammonia in 
dense molecular cores that find values as low as ours or even 
lower (e.g. Foster et al. [20091 Wootten [T995T l. 

4. Discussion 

4.1. Source multiplicity within OMC-2 FIR 4 



Several stu dies (J 0rgensen et al 120061 SOS, Grimier et al. 120091 
Kama et al. l2010l) suggest that OMC-2 FIR 4 contains three main 
components: 

- A large-scale cool envelope about 10000 AU (i.e. 0.05 pc) 
across. 

- An inner hot core, with kinetic temperatures above 100 K 
and around 400 AU in size. 

- A shock region produced by the interaction of an external 
bipolar outflow driven by the nearby OMC-2 FIR 3, which 
lies to the north-east of FIR 4. 

In addition, SOS report the presence of 1 1 cores seen at 3- 
mm inside this region, and speculate that such fragmentation has 
been externally triggered by the FIR 3 outflow. 

The emission of all the maps presented in Sect. |3] covers 
roughly 10000 AU (0.05 pc) and therefore traces the envelope 
component of OMC-2 FIR 4. In particular, this component ap- 
pears to be better traced by the NH3 maps, which are essentially 
single-peaked and display the most extended emission, similar 
in size and morphology to the structures detected in previous 
single-dish observations (e.g. Chini et al. 1 1 9971 the ASTE maps 
of Shimajiri et al. 200S). The temperature derived for the am- 
monia emission, ~ 30 K, is typical of protostellar (in our case 
protocluster) envelopes (e.g. Olmi et al.'2012'). 

At the resolution of our PdBI observations, despite being 
slightly better than that of the NMA observations, we cannot 
confirm the existence of the 1 1 cores claimed by SOS. Inspection 
of the different maps presented in this work allows us to clearly 
separate three regions, denoted by the identifiers main, west, and 
south (Sect. 13.2.1b . It appears, however, that our 2-mm contin- 
uum and some molecular line maps indeed hint at the presence 
of multiple smaller condensations within the three regions we 



detect, but these will need to be confirmed by higher angular 
resolution observations. 

At this point, it is worth mentioning that the hot core hypoth- 
esis proposed by Crimier et al. (2009) is based on a single source 
at the centre of the large-scale envelope. SOS and our new obser- 
vations show that the envelope component breaks down at scales 
of ~ 2000 AU, much larger than the presumed hot core size 
(with a radius of 440 AU). Only high-angular resolution imag- 
ing of typical hot core tracers (e.g. complex organic molecules) 
can observationally prove the existence of a hot core. 

As shown in Fig. |6] a velocity difference of about 6 km s"' 
can be seen between the south-west and the north-east extremes 
of the main source in the CH3OH and C^'*S(3-2) maps, with the 
north-eastern tip being the most red-shifted. SOS found the same 
behaviour in their single-dish methanol maps. Additional inter- 
ferometric observations of CO and SiO lines provided the same 
authors with morphological and dynamical evidence which led 
them to suggest that there is an interaction between the noith- 
eastern edge of OMC-2 FIR 4 and the molecular outflow driven 
by the nearby FIR 3. The red-shifted emission we detect from 
CH3OH and C^"^S(3-2) may also be explained by this scenario, 
i.e., it arises from shocked gas from the interaction between the 
FIR 3 outflow and the dense gas of FIR 4. It might alternatively 
be tracing the FIR 3 outflow itself, which cannot be detected 
much beyond the boundary of the FIR 4 clump due to the lim- 
ited sensitivity of the PdBI to more extended structures. More 
observations are necessary to confirm the external outflow inter- 
pretation, since with our data alone we are not able to conclude 
that this is indeed the case. Other possibilities exist that could 
cause a velocity gradient like the one we detect, such as rotation 
or internal outflowing/expanding motions. 

4.2. Nature of the sources 

Our observations allow us to distinguish three separate sources 
which may or may not present further substructure: 

1. The main source, at the nominal coordinates and systemic 
velocity ( 1 1 .4 km s" ' ) of OMC-2 FIR 4. It is the strongest 2- 
mm continuum and methanol source in the region. Its elon- 
gated shape suggests that it harbours two or more individual 
cores, although this hypothesis can only be confirmed with 
higher spatial resolution imaging. The red-shifted velocity 
of its northern portion might indicate an interaction with the 
red-shifted lobe of the external outflow driven by FIR 3. 

2. The west source, with a systemic velocity of 9-10 km s"'. It 
is weU traced by C^'*S(3-2) and DCN(2-1), and weaker than 
the main source emission. 

3. The south source, also with Vlsr =9-10 km s"', best traced 
by DCO+(2-l). It is the faintest methanol, DCN, and 2- 
mm continuum source, and has no significant emission in 
C'''*S(3-2). Its association with DCO^ indicates this region 
is likely colder than the other two, suggesting that it contains 
either a starless core or a protostar in an earlier evolution- 
ary stage or with a lower mass (e.g. J0rgensen et al. 12004) . 
Alternatively, since it does not display any emission peak in 
the 2-mm continuum image, it may not be a separated molec- 
ular core, but part of the cold external envelope of OMC- 
2 FIR 4. 

From the results of our LVG analysis on the CH3OH lines, 
we find that the entire OMC-2 FIR 4 clump is rather dense (« > 
2 X 10'' cm"^ at 15", i.e. 6000 AU). This suggests that the three 
regions we identify in this work, which are altogether smaller 
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Table 8. Results of the line fitting to the NH3 (1,1) and (2,2) lines and derived physical parameters 



Results of the hyperfine line fitting 


Line 


TsB 

(K) 


Vlsr 
(kms-') 


AV 
(kms-') 


''"main 




NH,(1,1) 
NH3(2,2) 


7.5 (2.7) 
4.9 (9.8) 


11.1(0.1) 
11.2(0.1) 


1.5(0.1) 
1.9(0.1) 


0.7 (0.2) 
0.1 (0.2) 




Physical parameters 


Tex 

(K) 


J rot 

(K) 


(K) 


A'nh, 
(10'*cm-2) 


Abu. (20 K)" 

(io-'») 


Abu. (50 K)" 
(lO-'O) 


18.2 


22.7 


29.6 


2.6 


1.4 


3.8 



M 



Abundance determined using the dust temperature indicated in brackets to derive the total gas column density (see Table[5j 



in size when compared to the whole emission considered in the 
CH3OH analysis, may be even denser. We note that the density 
we derive is in good agreement with SOS and with the predictions 
made by Grimier et al. (12009) , both of which report a value of 
10'' cm^l 



4.3. An ionised region in OMC-2 FIR 4 

Radio-continuum observations of the OMC2/3 filament were 
carried out with the VLA by Reipurth et al. (11999) at 3.6 cm. At 
the position of OMC-2 FIR 4, they detected an elongated source 
which the authors interpret as a radio jet. When comparing their 
map with our PdBI map, despite the different angular resolutions 
(i.e. ~ 4" for the PdBI maps, 8" for the VLA maps), we find 
that the position and shape of the radio source coincides fairly 
well within the pointing errors of the two maps, with the emis- 
sion distribution of our 2-mm continuum and methanol PdBI im- 
ages, which are both elongated in the same direction as the VLA 
source (see Fig.fTOt. 



While Reipuith et al. (119991 ) attribute the radio continuum 
emission to a jet, the interferometric observations performed to- 
wards this source so far (e.g. SOS, this work) have failed to detect 
molecular line emission at high velocities indicative of a jet or 
outflow emanating from this region. This leads us to consider a 
different possibility, namely that the radio emission arises from 
a very compact gas cavity ionised by a central young star, i.e. 
an Hii region, which is either surrounded by a shell of dust and 
molecular gas best traced by the 2-mm continuum and CH3OH 
emission, or else harbours two or more unresolved molecular 
dense cores within it. The Hii region scenario is a plausible ex- 
planation for the continuum radio emission, since IM protostars 
can also produce sufficiently intense Hii regions to be easily de- 
tected at centimetre wavelengths (e.g. Thompson , 1984) . 



Reipurth et al. CI 999) measured a flux of 0.64 mJy at 3.6 cm 
towards OMC-2 FIR 4. In the scenario of an optically thin Hii 
region, and assuming a single zero-age main sequence star is 
powering it, this implies a B3-B4 spectral type for the central 
star (using the tables of e.g. Panagia 119731 Thompson 119841 
Martins et al. 12005 D . This translates into a luminosity between 
700 and 1000 Lq and would indicate that this young star is the 
most massive source in the protocluster. It is also possible, given 
the elongated shape of the radio emission, that two (or more) in- 
dividuals are responsible for the ionisation of their surroundings. 
Only higher angular resolution and sensitivity images of the con- 
tinuum emission at different wavelengths will allow us to verify 
this possibility and find out whether the nature of the centimetre 
emission is an Hn region, as proposed here, or a radio jet. 



4.4. A final note on luminosities 

There is a certain confusion regarding the luminosity of OMC- 
2 FIR 4 in the literature. While the first estimates by Mezger et 
al. ( II99OI 1, and later by Crimier et al. (12009 ), report a value be- 
tween 400 L0 and 1000 Lq, a recent work by Adams et al. (12012) 
yields 50 Lq. In what follows we try to understand the reason of 
this discrepancy. 

The first point to mention is the complexity of the region. 
Indeed, the morphology and size of OMC-2 FIR 4 change when 
seen at different wavelengths: in the mid-IR it appears to be a 
compact source, while at wavelengths above ~ 100 ;um its emis- 
sion is more extended (~ 30" in diameter) and peaks around 
5" west of the mid-IR source. This said, we note that the lumi- 
nosity estimates published by the different authors refer to dif- 
ferent components. On the one hand, Adams et al. (2012) built 
the Spectral Energy Distribution of the source by measuring the 
continuum fluxes, at IR and sub-mm wavelengths, of the com- 
pact source they detect at 19 //m and 37 //m with SOFIA. In 
other words they adopted small aperture sizes, ranging between 
7" and 13" and centred at the mid-IR peak. On the other hand, 
Crimier et al. (2009) and Mezger et al. (11990) took into account 
a larger area, of 30" and 50" in diameter, respectively, which 
includes the whole envelope of OMC-2 FIR 4. In particular: 



- Crimier et al. (l2009l l used far-IR data from ISO and IRAS, 
whose coarse angular resolution likely led to measured 
fluxes that were overestimated, due to possible contamina- 
tion from the nearby protostar FIR 3, by at most a factor 2, 
since at such wavelengths the dominant emission source is 
FIR 4. 

- Adams et al. (120 1 21 ) obtained //e«c/zeZ-PACS maps of OMC- 
2 FIR 4 from which they measured the fluxes at 70 jum and 
160 /urn. They did so by applying a small aperture of 9.6" 
and 12.8", respectively, around the mid-IR peak coordinates, 
and a background annulus that actually covers part of the ex- 
tended emission of the region. In addition, they carried out 
APEX observations at 350 fim and 850 yum which lacked 
the sufficient sensitivity to detect any emission. They there- 
fore measured an upper limit of the flux considering only the 
emission within a beam size (7.3" and 19", respectively), 
again around the mid-IR coordinates, and therefore offset by 
~ 5" from the actual sub-mm peak position. As a result, their 
reported sub-mm fluxes are respectively 10 and 3 times lower 
than those detected and measured by Crimier et al. (2009) 
with the CSO and JCMT telescopes. As a result, considering 
these points, the authors may have underestimated the flux 
of OMC-2 FIR 4 by a factor of up to 10. 

In short, it appears that the discrepancy in luminosity esti- 
mates found in the literature is caused by the fact that different 
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Fig. 10. Left panel: VLA map of the 3.6 cm continuum emission in the OMC-2 FIR 4 region and its surroundings (taken from 
Reipurth et al. 1999). The synthesised beam is shown at the bottom-right comer. The source labelled VLA 12 coincides with the 
position of OMC-2 FIR 4. Right panel: Close-up view showing our PdBI continuum (grey scale) and CH3OH (black contours) 
maps. The synthesised beam is shown at the bottom-left corner. The red ellipse schematically represents the Act contour of the VLA 
image (left panel). The red cross marks the position of the cm emission peak. 



Table 9. Mass and luminosity estimates for OMC-2 FIR 4° 



Ref. 



L 


M 


Size 


(Lo) 


(Mo) 


(") 


350 


34 


50 


920 


30 


30 


50 


10 


8 



Mezgeretal. (1990} 
Crimieretal. (20093 
Adams etal. (2012) 



*"' Values re-scaled to a distance of 420 pc 



authors measure different components of OMC-2 FIR 4, with 
Adams et al (2012) concentrating on the compact mid-IR source 
and Crimier et al. (2009) on the entire envelope, which is sev- 
eral times larger. Table |9]summaiises the different values for the 
luminosity, mass and size of OMC-2 FIR 4 reported in the litera- 
ture. Based on the discussion above, the actual value depends on 
how one defines OMC-2 FIR 4 and its boundaries. If our inter- 
pretation of an Hii region is correct (see Sect. 14. 3) . the luminosity 
of the whole protocluster could be around several hundred solar 
luminosities. 



5. Conclusions 

We have presented new interferometric observations carried out 
towards the intermediate-mass protocluster OMC-2 FIR 4, in 
the Orion A complex. These include PdBI maps of the contin- 
uum emission at 143.4 GHz (2 mm) and the CH3OH(3,l,+0 
- 2,l,H-0), C^'*S(3-2), DCN(2-1), and DCO+(2-l) lines, as 
well as VLA maps of the NH3(1,1) and (2,2) inversion tran- 
sitions. Despite the coarse channel spacing of the observations 
(0.6 km s ' for the VLA spectra and 4 km s ' for the PdBI im- 
ages), we are able to extract important information about the 
morphology and kinematics of the region. Figure [TT] shows a 
cartoon of the OMC-2 FIR 4/FIR 3 region, summarising all the 
elements present according to our interpretation of the data pre- 
sented in this work. Our findings are summarised as follows: 



1 . Our 2-mm continuum maps reveal two main components and 
hints of secondary cores or clumps surrounding it. The exis- 
tence of the 1 1 cores discovered by SOS cannot be confirmed 
by our observations, despite our higher angular resolution. 

2. Each PdBI line traces the region differently. By compar- 
ing their emission, we are able to distinguish three sepa- 
rate sources of one or several solar masses each, which we 
have named main, west and south. The main and west ones 
are warmer. The south source displays a stronger line emis- 
sion in DCO^(2 - 1), suggesting a colder temperature and 
therefore an earlier evolutionary state or lower mass core. 
Alternatively, it may be associatied with part of the colder 
envelope of OMC-2 FIR 4. As for the ammonia maps, their 
emission is more extended and essentially single-peaked, 
and therefore traces the envelope of the protocluster. An 
LVG modelling of the CH3OH lines we observe allows us 
to conclude that the whole protocluster is very dense, with 
n > 2 X 10* cm""*, in agreement with previous theoretical 
predictions on the structure of the envelope of OMC-2 FIR 4. 

3. The main source as seen by the 2-mm continuum and 
CH3OH emission coincides, both in position (within the 
pointing eiTors) and in shape, with the radio-continuum 
emission observed at 3.6 cm with the VLA (Reipurth et 
al. 119991) . In absence of clear interferometric evidence of a 
molecular outflow driven by FIR 4, we interpret this as an 
Hi: region powered either by a B3-B4 type young star, or by 
two or more lower-mass stellar objects. This region is likely 
suiTounded by a shell of dust and molecular gas (seen by the 
millimetre continuum and CH3OH maps). This scenario re- 
mains to be verified by higher angular resolution and multi- 
wavelength observations. 

Due to its relative proximity to the Sun, OMC-2 FIR 4 is an 
ideal laboratory to explore intermediate- and low-mass star for- 
mation in the context of a whole protocluster. The evidence of 
core multiplicity in this region, coupled with the insufficient spa- 
tial resolution of our observations to clearly separate them, high- 
lights the need of higher spectral and angular resolution imag- 
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Fig. 11. Schematic cartoon of OMC-2 FIR 4 and the nearby 
OMC-2 FIR 3, summarising all the physical elements in the 
region that we propose according to the data presented in this 
work (see Sect. |4|l. The white ellipse inside the main source de- 
picts the ionised region powered by a B4 young star. The colours 
of the main, west, and south sources represent their systemic 
velocities with respect to the nominal value for OMC-2 FIR 4 
(Vlsr = 1 1.4 km s"'), marked in green. The small red and blue 
circles within the three sources represent the possibility that they 
harbour smaller unresolved molecular condensations, with the 
colour blue denoting a colder temperature than red. 
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ing to disentangle the different components and extract more 
information about the kinematics and the excitation conditions 
in each of them. This information will be of crucial importance 
to interpret the Herschel-HIFI data obtained towards this source 
within the CHESS Key Programme. 
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